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The  goal  of  this  research  project  Is  to  study  the  use  of  the 
polarization  properties  of  the  signal  scattered  from  a  radar  target  for 
the  purpose  of  radar  target  classification  or  identification.  The 
research  Involves  the  application  of  a  large  data  base  derived  from 
scale  model  measurements  using  the  Ohio  State  University  Compact  Radar 
Measurement  Range.  The  data  base  has  been  combined  with  a  set  of 
polarization  state  transformation  algorithms  to  permit  the  testing  of 


NCTR  algorithms  as  a  function  of  polarization.  Also  studied  is  the 
application  of  complex  natural  resonance  techniques  to  the  scattering 
data  for  various  polarization  states. 

It  is  shown  that  co-polarized  fully  coherent  radar  systems  with 
multi  frequency  capability  can  perform  very  effectly  as  NCTR  systems. 
Signal -to-noise-ratlos  from  0  to  +10  dB  in  the  final  estimate  of  the  RCS 
and  relative  signal  phase  for  the  various  targets  is  shown  to  be 
sufficient  for  mlsclassification  percentages  less  than  5  percent.  It  is 
also  shown  that  natural  resonance  data  can  be  meaningfully  utilized  when 
applied  to  resonance  region  RCS  data. 
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INTRODUCTION 

*■»/ 

This  report  describes  the  results  of  the  first  year  of  a  research 
project  performed  by  the  Ohio  State  University  ElectroSclence  Laboratory 
(OSU/ESL)  for  the  Naval  Weapons  Center  (NWC).  The  goal  of  this  project 
Is  to  explore  the  use  of  the  polarization  properties  of  the  signal 
scattered  from  a  radar  target  for  the  purpose  of  radar  target 
Identification.  Various  radar  target  Identification  algorithms  were 
applied  to  the  case  of  a  full  polarization  coherent  radar  system,  and 
were  tested  using  a  specific  data  base  and  noise  model. 

The  data  base  used  to  test  the  performance  of  the  radar  target 
Identification  algorithms  developed  here  is  a  unique  set  of  measurements 
made  on  scale  models  of  aircraft.  Measurements  were  made  using  the 
OSU/ESL  Compact  Radar  Measurement  Range^tll.  The  range  was  operated  In 
a  broad-band  (1-12  GHZ)  mode  and  the  full  polarization  matrix  was 
measured.  Calibrated  values  (amplitude  and  phase)  of  the  RCS  for  the 
three  polarization  states  were  thus  available.  The  polarization  states 
are  listed  below. 
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INTRINSIC  POLARIZATION  STATES  OF  THE  RCS  MEASUREMENT  SYSTEM 

1.  VV... transmit  vertical  polarization;  receive  vertical 

polarization 

2.  HH.. .transmit  horizontal  polarization;  receive  horizontal 

polarization 

3.  HV.. .transmit  horizontal  polarization;  receive  vertical 

polarization 

The  models  used  .wrc  plasti  scale  models  of  modern  commercial 
aircraft  which  were  si  er  coated  using  a  dipping/plating  process.  The 
scale  factor  was  approximately  200,  which  put  the  full  scale  frequency 
range  In  the  5  MHZ  to  60  MHZ  frequency  band.  This  frequency/scale- 
factor  relationship  results  In  RCS  values  for  the  resonance  region  of 
the  various  aircraft.  Aspect  angles  from  0  degrees  (nose-on)  to  180 
degrees  (tall -on)  were  used. 

One  of  the  first  tasks  was  to  develop  the  polarization  state 
transformations  and  Incorporate  them  Into  a  computer  data  base 
manipulation  program  so  that  the  polarization  state  could  be  derived 
from  the  original  measurements  and  the  resultant  entered  Into  radar 
target  Identification  algorithms.  This  program,  called  RSSE,  Is 
described  In  ESL  technical  report  716559-1;  which  was  prepared  under 

this  program  to  document  the  algorithms  developed  during  this  research 
[2],  During  this  phase  of  the  research,  a  computer  program  which  would 

display  the  polarization  state  of  the  signal  scattered  from  the  radar 

A  A 

target  on  the  well-known  Poincare  sphere  w*s  developed.  The  Poincare 
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projections  for  the  various  aircraft  at  various  aspect  angles  and  for 
various  Illumination  polarization  states  proved  to  be  very  useful  for 
gaining  Insight  Into  the  utilization  of  such  a  mapping  for  the 
development  of  other  radar  Identification  algorithms.  Examples  are 
given  In  Chapter  3. 

The  application  of  resonance-region  radar  target  Identification 
algorithms  to  the  polarization-derived  features  Is  also  demonstrated  In 
Chapter  3.  Tests  of  the  performance  of  the  various  features  In  the 
radar  target  Identification  environment  are  shown  and  the  results 
statistically  compared.  These  comparisons  of  the  various  algorithms 
which  utilize  the  various  polarization  dependent  feature  sets  are  used 
to  summarize  the  results  of  this  phase  of  the  research  In  this  area. 

The  other  area  which  was  studied  Is  the  application  of  complex 
natural  resonance  techniques  to  polarization  states  and  frequency  bands 
which  were  not  available  In  the  previous  work.  This  research  Involves 
the  extraction  of  the  complex  poles  of  the  RCS  signatures  of  the  radar 
targets  described  In  Chapter  2  of  this  report.  The  extracted  poles  are 
then  used  In  specific  radar  target  Identification  algorithms  [3]. 
Statistics  of  the  performance  of  these  algorithms  are  given  In 
Chapter  3. 
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CHAPTER  2 


DATA  BASE  AND  HANDLING 

2.1  DATA-BASE  DESCRIPTION 

A  critical  prerequisite  for  performing  the  target  Identification 
research  discussed  below  Is  the  generation  and  malntanence  of  a 
well -organized,  easily  accessible  data  base  containing  the  precise 
wlde-band  full -polarization  scattering  matrix  measurements  of  a  number 
of  aircraft  targets  of  Interest.  Due  In  large  part  to  support  from  the 
Office  of  Naval  Research  (under  contract  N00014-82-K-0037) ,  such  a  data 
base  has  been  generated  at  the  ElectroSclence  Laboratory  using  the 
compact  radar  range  at  ESL  to  measure  the  far-fleld  radar  backscatter 
returns  from  scale-models  of  five  commercial  aircraft  Including:  the 
Concorde,  Douglas  DC-10,  and  the  Boeing  707,  727,  and  747. 

The  cross  section  measurements  made  on  these  scale-model  aircraft 
were  typically  performed  in  the  1-12  GHz  frequency  band  which,  for  a 
scale  factor  of  200:1  (say)  corresponds  to  full-scale  measurements  from 
5-60  MHz.  These  measurements  were  performed  for  each  of  the  aircraft 
listed  above  for  a  number  of  azimuth  angles  of  Interest,  using  both 
linear  co-pole  configurations  (HH  and  VV)  as  well  as  the  linear 
cross-pole  (HV)  configurations  of  the  compact  range.  Listings  of  the 
data-base  entries  for  each  aircraft,  Including  the  data-base  file 
locations,  the  scale  factor  of  the  model  employed,  and  the  range  of 


available  measurements  (In  GHz)  for  each  polarization  and  azimuth 
(aspect)  angle  appear  in  Tables  1-5  [2],  In  these  tables,  the  entry 
"6-12“  (for  example)  Indicates  that  only  measurements  in  the  6-12  GHz 
range  are  avlalable  for  this  aircraft  at  the  particular  polarization  and 
azimuth  angle.  The  entry  "null"  Implies  that  no  data  are  available  for 
this  target  at  this  azimuth  and  polarization. 

2.2  DATA-BASE  HANDLING 

Three  software  packages  for  data-base  manipulation  and  processing 
were  developed  under  support  or  partial  support  from  this  program;  each 
performing  a  set  of  tasks  of  Interest.  The  first  package  performs  tasks 
related  to  data-base  management.  Including  generation,  mapping,  and  any 
necessary  modifications.  The  second  performs  polarization  mappings  and 
Fourier  analyses  of  original  and  transformed  parameters.  The  third 
software  package  developed  under  the  sponsorship  of  this  program  allows 
the  Investigation  of  the  performance  of  various  radar  target 
Identification  systems  based  on  different  types  and  ranges  of  features. 

A  complete  description  of  the  use  and  capabilities  of  these 
computer  subroutines  appears  In  the  earlier  report  from  this  project 
[2],  Much  of  the  discussion  and  data  presented  below  are  based  on 
results  obtained  using  one  or  more  of  these  software  packages. 


TABLE  2.1 


DATA  BASE  ENTRIES  FOR  BOEING  707 


NAPPING  DATABASE  i 


TITLE  i  707  AIRPLANE 


CAL  I  MATED  FILES  FROM  DISK  PACK  0NR7  DIRECTORY  707 
DBM  HEADER  FILE  DLAOtEONR7.7O70707CAL.DBM 
SCALE  FACTOR  ■  190.0  ELEVATION  ANGLE  •  0  DECREES 


FREQUENCY  FORMATTED  DATA  BASE  (  > 


•  HH  • 

•  HV  • 

•  VV  • 

0 

1-12 

1-12 

1-12 

10 

1-12 

NULL 

1-12 

15 

1-12 

1-12 

1-12 

20 

1-12 

1-12 

1-12 

25 

6-12 

NULL 

6-12 

50 

1-12 

1-12 

1-6.3 

35 

6-12 

NULL 

6-12 

40 

1-12 

1-12 

1-12 

45 

1-12 

1-12 

1-12 

50 

1-12 

1-12 

1-12 

55 

6-12 

NULL 

6-12 

60 

1-12 

1-12 

1-12 

65 

6-12 

NULL 

6-12 

70 

1-12 

1-12 

1-12 

75 

1-12 

1-12 

1-12 

ac 

1-12 

1-12 

1-12 

85 

6-12 

NULL 

6-12 

90 

1-12 

1-12 

1-12 

95 

6-12 

HILL 

6-12 

100 

1-12 

1-12 

1-12 

105 

1-12 

1-12 

1-12 

110 

1-12 

1-12 

1-12 

115 

6-12 

HILL 

6-12 

120 

1-12 

1-12 

1-12 

125 

6-12 

NULL 

6-12 

130 

1-12 

1-12 

1-12 

135 

1-12 

1-12 

1-12 

140 

1-12 

1-12 

1-12 

145 

6-12 

NULL 

6-12 

150 

1-12 

1-12 

1-12 

153 

6-12 

NULL 

6-12 

160 

1-12 

1-12 

1-12 

165 

1-12 

1-12 

1-12 

170 

1-12 

1-12 

1-12 

173 

6-12 

NULL 

6-12 

180 

1-12 

1-12 

1-12 

270 

6-12 

NULL 

1-12 

'  <*.  .  /.V.  iCyVw  <Cg  V.V.Vy  V.V  j  < 


TABLE  2.2 


DATA  BASE  ENTRIES  FOR  BOEING  727 


RAPPING  DATABASE  I 
727 


TITLE  i  72?  AIRH.MC  f  INCH  727 

CA1 1  BRA  TED  FILES  FR»  DISK  FAOt  0NR1  DIRECTORY  727 
DBM  ICADCR  FILE  OLAOs[ONR1  .727 J7 27 CAL. DBM 
SCALE  FACTOR  >  200.0  ELEVATION  ANGLE  -  0  DEGREES 


FREQUENCY  FORMATTED  DATA  BASE  (  GHz  ) 


ASPECT  (  0*g 


POLARIZATION 


•  MM  • 

•KY¬ 

0 

1-12 

I-12 

10 

1-12 

15 

1-12 

1-12 

20 

1-12 

1-12 

30 

1-12 

1-12 

40 

1-12 

1-12 

45 

1-12 

1-12 

50 

1-12 

1-12 

60 

1-12 

1-12 

70 

1-12 

1-12 

75 

1-12 

1-12 

eo 

1-12 

1-12 

GO 

1-12 

1-12 

100 

1-12 

1-12 

105 

1-12 

1-12 

110 

1-12 

1-12 

120 

1-12 

1-12 

130 

1-12 

1-12 

135 

1-12 

1-12 

140 

1-12 

1-12 

150 

1-12 

1-12 

160 

1-12 

1-12 

165 

1-12 

1-12 

170 

1-12 

1-12 

160 

1-12 

1-12 

270 

1-12 

1-12 

•  yy  i 
1-12 
1-12 
1-12 
1-12 
1-12 
1-12 
1-12 
1-12 
1-12 
1-12 
1-12 
1-12 
1-12 
1-12 
1-12 
1-12 
1-12 
1-12 
1-12 
1-12 
1-12 
1-12 
1-12 
1-12 
1-12 
1-12 


TABLE  2.3 

DATA  BASE  ENTRIES  FOR  BOEING  747 


NAPPING  DATABASE  • 
747 


TITLE  i  747  AIRPLWC  LARGE  747 

CALIBRATED  FILES  FROM  DISK  PACK  0NR4  DIRECTORS  747 
DBM  ICADER  FILE  DLA0t[CNR4.7473747OAL.DGN 
SCALE  FACTOR  -  200.0  ELEVATION  ANGLE  ■  0  DEGREES 


FREQUENCY  FORMATTED  DATA  BASE  (  GHz  ) 


ASPECT  (  0*0  1 


0 

•  HH  • 
1.3-12 

10 

1.3-12 

15 

6-12 

20 

1.3-12 

25 

6-12 

30 

1.3-12 

35 

6-12 

40 

1.3-12 

45 

1.5-12 

30 

1.3-12 

33 

6-12 

60 

1.5-12 

65 

6-12 

70 

1.3-12 

75 

6-12 

60 

1.3-12 

GO 

1.3-12 

95 

6-12 

100 

1.3-12 

105 

6-12 

110 

1  3-12 

113 

6-12 

120 

1.5-12 

123 

6-12 

130 

1.3-12 

133 

6-12 

140 

1.3-12 

145 

6-12 

130 

1.3-12 

153 

6-12 

160 

1.3-12 

165 

6-12 

170 

1.3-12 

iec 

1.3-12 

POLARIZATION 


•  HV  • 

•  VV  • 

1-12 

1-12 

1-12 

NULL 

1-12 

1-12 

1-12 

NULL 

NULL 

NULL 

1-12 

1-12 

NULL 

NULL 

1-12 

NULL 

1-12 

1-12 

1-12 

NULL 

NULL 

MILL 

1-12 

1-12 

NULL 

MILL 

1-12 

NULL 

1-12 

1-12 

1-12 

NULL 

1-12 

1-12 

NULL 

NULL 

1-12 

MILL 

1-12 

1-12 

1-12 

NULL 

NULL 

MILL 

1-12 

1-12 

MILL 

NULL 

1-12 

NULL 

1-12 

1-12 

1-12 

NULL 

MILL 

NULL 

1-12 

1-12 

MILL 

NULL 

1-12 

NULL 

1-12 

1-12 

1-12 

NULL 

1-12 

1-12 
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TARLE  2.4 


DATA  BASE  ENTRIES  FOR  THE  CONCORDE 


NAPPING  DATABASE  I 
CON 


TITLE  t  CONOOM)  A1RHAIC  URGE  CONCORD 

CALIBRATED  FILES  FROM  DISK  RACK  ONR1  DIRECTORY  CON 
DBM  HEADER  FILE  DLAO>[ONR1.CON3CONQAL.OBM 
SCALE  FACTOR  -  130.0  ELEVATION  ANGLE  •  0  DEGREES 


FREQUENCY  FORMATTED  DATA  BASE  (  GHz  ) 

ASPECT  (  Dag  )  POLARIZATION 


•  MH  • 

■  Hi  ■ 

•  TV  • 

0 

1-12 

1-12 

1-12 

10 

1-12 

1-12 

1-12 

15 

1-12 

NULL 

1-12 

20 

1-12 

1-12 

1-12 

25 

8-12 

NULL 

6-12 

50 

1-12 

1-12 

1-12 

55 

6-12 

NULL 

6-12 

40 

1-12 

1-12 

1-12 

45 

1-12 

1-12 

1-12 

50 

1-12 

1-12 

1-12 

55 

6-12 

NULL 

6-1* 

SO 

1-12 

1-12 

1-12 

S3 

6-12 

NULL 

6-12 

70 

1-12 

1-12 

1-12 

75 

1-12 

1-12 

1-12 

BO 

1-12 

1-12 

1-12 

85 

6-12 

NULL 

6-12 

90 

1-12 

1-12 

1-12 

95 

6-12 

MILL 

6-12 

100 

1-12 

1-12 

1-12 

105 

1-12 

1-12 

1-12 

110 

1-12 

1-12 

1-12 

115 

6-12 

MILL 

6-12 

120 

1-12 

1-12 

1-12 

123 

6-12 

NULL 

6-12 

130 

1-12 

1-12 

1-12 

133 

1-12 

1-12 

1-12 

140 

1-12 

1-12 

1-12 

145 

6-12 

NULL 

6-12 

150 

1-12 

1-12 

1-12 

155 

6-12 

NULL 

6-12 

160 

1-12 

1-12 

1-12 

165 

1-12 

1-12 

1-12 

170 

1-12 

1-12 

1-12 

175 

6-12 

NULL 

6-12 

180 

1-12 

1-12 

1-12 
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TABLE  2.5 

DATA  BASE  ENTRIES  FOR  DCIO 


MAPPING  DATABASE  t 
DCIO 


TITIE  i  DCIO  AIRPLAIC  ELEVEN  INCH  DCIO 

CALIBRATED  FILES  FROM  DISK  PACK  0NR7  DIRECTORY  DCIO 
DBM  HEADER  FILE  DLAOi[0NR7.DC10]DC1OCAL.DBM 
SCALE  FACTOR  -  200.0  ELEVATION  ANGLE  >  0  DESREES 


FREQUENCY  FORMATTED  DATA  BASE  (  GHz  I 

ASPECT  (  D*S  >  POLARIZATION 


•  MM  • 

»  MV  ■ 

•  VV  ■ 

to 

1-12 

1-12 

1-12 

15 

1-12 

1-12 

1-12 

20 

1-12 

1-12 

1-12 

30 

1-12 

1-12 

1-12 

40 

,  1-1? 

1-12 

1-12 

43 

1-12 

1-12 

1-12 

SO 

1-12 

1-12 

1-12 

60 

1-i* 

1-12 

1-12 

70 

1-12 

1-12 

1-12 

75 

1-12 

1-12 

1-12 

BO 

1-12 

1-12 

1-12 

00 

1-12 

1-12 

1-12 

100 

1-12 

1-12 

1-12 

105 

1-12 

1-12 

1-12 

110 

1-12 

1-12 

1-12 

120 

1-12 

NULL 

1-12 

130 

1-12 

NULL 

1-12 

135 

1-12 

NULL 

1-12 

140 

1-12 

1-12 

1-12 

150 

1-12 

1-12 

1-12 

160 

1-12 

1-12 

1-12 

163 

1-12 

1-12 

1-12 

170 

1-12 

1-12 

1-12 

180 

1-12 

1-12 

1-12 

270 

1-12 

1-12 

1-12 
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FULL  SCATTERING  MATRIX  TARGET  ID  RESEARCH 

3.1  INTRODUCTION 

Several  Investigators  have  suggested  that  data  containing  the 
broad-band  measurements  of  the  full  polarization  scattering  matrix  for  a 
set  of  targets  would  contain  a  great  deal  of  Information  useful  for 
target  Identification  algorithms  [41;  however,  until  recently,  precise 
measurements  for  all  but  canonical  targets  were  not  available  and  target 
Identification  techniques  based  on  polarization-intensive  Information 
could  not  be  evaluated  for  practical  radar  targets. 

The  primary  focus  of  the  research  described  In  this  report,  and  the 
motivation  for  the  generation  of  the  data  base  and  handling  software 
discussed  above,  was  the  investigation  of  the  performance  of  radar 
target  Identification  algorithms  applied  to  the  full  polarization 
scattering  matrix  data  for  a  set  of  five  common  aircraft.  As  a  part  of 
this  study,  the  effects  of  additive  white  Gaussian  noise  on  the 
performance  of  systems  employing  various  polarizations  and  combinations 
of  polarizations  was  Investigated.  Of  particular  Interest  were  the 
Initial  comparisons  of  the  performance  of  a  nearest-neighbor  algorithm 
using  linear  (HH,  VV,  and  HV),  circular  (RR,  LL,  and  RL),  and 
axial-ratio  polarization  coefficients  as  features. 
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A  secondary  goal  of  this  study  was  the  assessment  of  the  value  of 
precise  phase  Information  In  this  particular  additive  noise  environment. 
This  was  accomplished  by  performing  simulation  studies  using  both 
"coherent"  (complex)  and  "non-coheren*"  (magnitude)  distance  metrics  as 

t 

part  of  the  nearest-neighbor  algorithm.  In  each  case,  the  performance 
parameter  of  Interest  was  the  average  percentage  of  mlsclasslflcatlon  as 
a  function  of  the  average  power  level  of  the  noise. 

In  the  following  sections  of  this  chapter,  several  of  the  key 
issues  In  polarization  scattering  matrix  based  target  Identification  are 
discussed  In  more  detail  and  the  results  obtained  as  part  of  this  effort 
are  presented  and  compared. 

3.2  POLARIZATION  STATE  TRANSFORMATIONS 

The  data  base  contains  the  absolute  amplitude  and  phase  of  the  RCS 
of  five  different  aircraft  over  a  broad  band  of  frequencies. 

Measurements  were  made  using  vertical  and  horizontal  polarization  states 
( VV-HH-VH ) .  It  Is  thus  possible  to  transform  from  the  linear 
polarization  states  to  circular  (or  any  elliptical  polarization  state). 
Listed  below  are  several  reasons  for  considering  this  type  of 
transformation  of  the  data  prior  to  applying  a  target  Identification 
scheme. 

•  They  are  linear  transformations;  so  the  characteristics  of  the 
additive  noise  process  will  remain  unchanged. 
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•  Aircraft  and  missiles  In  particular  tend  to  be  symmetric 
bodies.  Considering  the  head-on  or  tall -on  aspect  angle,  It 
can  be  seen  that  a  change  In  the  roll  angle  of  the  aircraft 
will  change  the  vertical  and  horizontal  scattering  components 
considerably,  but  will  not  change  the  circular  scattering 
components  except  In  the  phase, 

•  The  symmetry  of  aircraft  and  missiles  also  tends  to  transform 
the  scattering  matrix  [T]£X2  t0  4  unitary  matrix  If  the 
components  are  circularly  polarized.  Observe  that  In  the  nose 
or  tall  aspect  zone,  the  symmetry  of  the  aircraft  or  missile 
tends  Lo  make  T 12  and  T21  equal  to  zero  (as  Is  also  the  case 
for  linear  components),  and  that  It  has  the  additional  property 
of  making  Tjj  equal  to  T22. 

•  Even  In  the  non-symnetrlc  aspect  zones,  aircraft  or  missile 
maneouvers  will  have  less  Influence  on  the  parameters  of  the 
circular  polarization  RCS  components  because  components  of 
rotation  about  the  line  of  sight  axis  will  not  enter  Into  the 
response  variation. 

With  these  considerations  In  mind,  a  polarization  state 
transformation  algorithm  was  implemented  and  comparisons  of  the  target 
Identification  algorithm  performance  were  made. 


3.2.1  Linear  To  Circular  Polarization  Transformation 
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In  general,  a  radar  system  Including  the  radar  target  can  be 
considered  a  two-port  network  problem.  The  relationships  can  be 
expressed  as 
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(3.2) 


Where  Sr  is  the  received  signal,  Sj  Is  the  transmitted  signal, 
and  Tjj  expresses  the  relationship  between  them.  In  our  case,  T^j  Is 
normalized  to  give  the  relationship  in  units  of  meters  (l.e.,  the  square 
root  of  the  RCS  In  units  of  square  meters).  Since  reciprocity  holds 
here,  T^  equals  T^. 

In  these  equations,  the  subscripts  1  or  2  Indicates  orthogonal 
polarization  states.  The  relationships  can  be  expressed  In  terms  of  any 
two  orthogonal  polarization  states.  There  exist  transformations  between 
the  various  orthogonal  polarization  states,  and  we  will  use  the 
transformations  from  linear  states  to  circular  states  here.  They  are 
given  as:  1 


\ 

IPR  =  2  (lw  “  IhH^  +  J'lVH 


(3.3) 
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(3.4) 
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III  "  2  <Ivv  "  Ihh^  ”  jIvh 
1 

IlR  *  IrL  ■  2  <Iw  +  W  (3.5) 

In  this  case,  the  subscripts  are  defined  as; 

*  R  «  right  hand  circular 

*  L  ■  left  hand  circular 

*  V  ■  vertical 

*  H  *  horizontal 

Where  Tjj  are  given  In  units  of  area  (RCS)  and  phase  with  respect  to  a 
reference  point  near  the  center  of  the  target. 

These  transformations  have  been  Implemented  on  the  ElectroScience 
Laboratory  VAX-11/7R0  computer  to  be  compatable  with  the  aircraft  data 
base.  The  transformed  results  are  then  used  as  Input  to  the  target 
Identification  algorithms.  This  means  that  the  Input  to  the  target 
Identification  algorithms  can  be  any  of  the  feature  sets  described 
below.  (The  program  Is  documented  In  detail  In  report  716559-1  [21.) 

Features  available  to  the  target  in  algorithm 

1.  Radar  Cross  Section  (calibrated  amplitude  and  phase) 

*  transmit  vertical  or  horizontal  and  receive  vertical  or 
horizontal  (VV  -  HH  -  VH  or  HV) 

*  transmit  right  or  left  hand  circular  and  receive  right  or 
left  hand  circular  (RR  -  LL  -  RL  or  LR) 
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2.  transmit  right  or  left  hand  circular  and  measure  the  axial 
ratio  of  the  received  polarization  ellipse  (RAX  -  LAX) 

Note  the  utilization  of  the  received  signal  polarization  axial 
r*at1o.  This  Is  computed  from  the  orthogonal  components  of  the  received 
signal.  If  the  received  signal  were  displayed  as  a  signal  voltage  In 
the  plane  perpendicular  to  the  direction  of  propagation,  the  signal 
would  trace  out  the  well-known  polarization  ellipse.  The  ratio  of  the 
major  to  minor  axis  of  this  ellipse  Is  the  axial  ratio  referred  to 
here. 

Note  that  the  polarization  nulls  are  not  used  as  a  target 
Identification  feature.  These  parameters  of  the  radar  scattering  from  a 
target  have  the  advantage  that  they  are  radar-independent  features  of 
the  target  (they  are  dependent  on  the  pitch,  roll  and  yaw  of  the 
target).  They  may  be  studied  at  some  future  date,  but  In  terms  of  their 
Information  content  and  usefulness  for  target  Identification,  there  Is 
no  Indication  at  this  time  that  they  will  contain  more  features  relevant 
to  target  Identification  than  the  direct  measurements  discussed  here. 

Examples  of  the  use  of  the  above  features  In  the  Identification  of 
the  aircraft  targets  In  the  data  base  will  be  shown  In  the  next 
section. 

3.2.2  Poincare  Transformations 

It  Is  possible  to  display  the  measurements  made  on  the  aircraft 
data  base  on  the  well-known  Poincare  sphere.  In  this  case,  the 
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polarization  state  of  the  received  signal  for  various  transmitted  signal 
polarizations  Is  presented  as  a  point  plotted  on  the  surface  of  a 

0 

sphere.  An  example  of  the  Poincare  sphere  Is  shown  In  Figure  3.1.  This 
figure  shows  the  relationship  between  a,  B,  y  and  $,  where  these 
parameters  are  listed  below. 

0 

Poincare  parameters 

*  a  »  arctangent  of  ellipse  axial  ratio 

*  B  *  orientation  of  ellipse  major  axis 

*  y  •  arctangent  of  the  ratio  of  the  vertical  and 
horizontal  components  of  the  signal  current  magnitudes 

*  ♦  *  electrical  phase  angle  between  the  vertical  and 
horizontal  signal  components 


An  example  where  the  scattering  from  the  Boeing  747  model  has  been 
transformed  to  show  the  variation  of  the  tilt  angle,  the  axial  ratio  and 
the  Intensity  of  the  scattered  signal  for  the  case  of  left-hand  circular 
Illumination  Is  shown  in  Figure  3.2.  Note  that  these  parameters  are 
shown  as  a  function  of  frequency  over  the  measurement  band  from  1  to  12 
GHz.  Note  that  In  this  case,  the  target  Is  viewed  nose-on.  This  means 
that  the  aircraft  Is  a  symmetric  radar  target,  and  If  the  target  had 
been  Illuminated  with  a  left-hand  circular  transmitted  signal,  the 
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Figure  3.1.  Polar  projection  of  Poincare  sphere 


Intensity  would  have  been  the  same,  but  the  axial  ratio  and  the  tilt 
angle  would  have  been  reversed.  (These  results  were  confirmed  with  the 
actual  measured  data.)  The  Impact  of  this  observation  Is  that  the 
number  of  Independent  variables  Is  significantly  reduced  In  this  feature 

t 

space.  There  Is  no  reason  to  transmit  the  orthogonal  circular 
polarization  state  for  the  nose-on  (or  tall-cn)  aspect  angle  because  the 
resulting  data  will  not  contain  any  new  Information.  If  the  roll  angle 
of  the  target  were  to  change,  only  the  tilt  angle  parameter  would 
change^  Furthermore,  this  parameter  would  change  In  a  predictable  way 
(It  would  roll  with  the  target).  From  the  point  of  view  of  *arget 
Identification,  this  means  that  the  Intensity  and  the  axial  ratio  are 
quite  stable  as  the  aircraft  rolls.  From  the  point  of  view  of  radar 
signature  analysis,  this  means  that  the  roll  angle  of  the  target  Is 
contained  In  the  tilt  angle  parameter.  This  may  be  a  useful  parameter 
for  predicting  the  future  position  of  the  aircraft  (or  missile)  since 
the  aircraft  Is  required  to  roll  in  order  to  Initiate  a  turn. 

A  Poincare  plot  of  the  parameteres  shown  In  Figure  3.2  Is  shown  in 
Figure  3.3.  Note  that  the  amplitude  Is  Indicated  by  the  size  of  the 
symbol  plotted  at  each  point.  The  data  proceeds  from  point  to  point  as 
a  function  of  frequency.  Once  again,  the  starting  frequency  was  1  GHz 
(marked  with  an  “S"  in  the  figure),  the  Increment  Is  50  MHz,  and  the 
final  frequency  Is  12  GHz  (marked  with  an  "F”  on  the  figure).  (The 
middle  frequency  (5.5  GHz)  Is  marked  with  an  "M1'  on  the  figure.) 

Note  that  from  the  perspective  of  a  human  being,  this  Is  a  messy 
plot.  Certain  topological  patterns  begin  to  emerge,  however,  as  the 
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plot  Is  studied.  Note  for  example  that  the  larger  amplitude  data  falls 
Into  three  clusters;  one  near  linear  horizontal,  one  near  linear 
vertical  and  one  near  linear  45  degrees  tilt  angle.  It  can  also  be  seen 
that  the  most  common  variation  as  a  function  of  frequency  Is  a  variation 
in  tilt  angle.  Automatic  algorithms  for  the  extractions  of  such 
parameters  and  their  utilization  In  target  Identification  schemes  will 
be  studied  In  the  near  future. 

The  key  to  the  utilization  of  such  features  for  target 
Identlflcaton  lies  In  the  stability  of  these  features  for  a  given 
aircraft  type  and  the  variation  of  these  features  between  different 
aircraft  types.  A  set  of  examples  where  this  property  can  be  evaluated 
Is  shown  In  Figures  3.4  through  3.8.  In  these  figures,  the  data  from 
five  aircraft  are  shown  plotted  as  In  Figure  3.3  but  with  a  much  smaller 
frequency  extent.  The  frequencies  were  chosen  so  as  to  cover  the  same 
full-scale  band  (the  8  to  16  MHz  band)  for  the  various  aircraft.  Note 
that  this  frequency  band  covers  a  bandwidth  of  8  MHz,  which  implys  a 
range  resolution  of  approximately  19  meters  (61  feet),  therefore,  we  are 
clearly  not  using  a  bandwidth  large  enough  to  resolve  individual 
scattering  centers.  The  target  Is  behaving  much  like  a  lumped  Impedance 
transformer  (1e:  It  generates  an  amplitude  and  phase  response  as  a 
function  of  frequency,  but  such  response  does  not  resolve  Into  discrete 
scattering  centers.) 

Note  that  the  patterns  seen  In  the  data  are  simpler  for  hunnn 
viewing  but  that  a  computer  may  be  just  as  happy  with  the  larger 

bandwidth  data  shown  In  Figure  3.3.  Also  note  that  the  data  tend  to 
cluster  the  high  amplitude  components  near  the  linear  polarization  state 
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Figure  3.2.  Polarization  parameters  for  Boeing  747  as  a  function  of 
frequency. 
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Polar  projection  of  Poincare  sphere  representation  of 
scattering  f;om  Boeing  747. 


for  all  of  the  models.  An  occasional  "looping"  pattern  can  also  be  seen 

In  some  of  the  data  (the  Concorde  and  the  727  for  example).  We  have 

found  that  this  Is  usually  due  to  Interference  between  multiple 

scattering  centers  on  a  target. 

« 

Of  particular  Interest  Is  the  practicality  of  using  these  types  of 
presentations  to  understand  the  behavior  of  the  polarization  state 
parameters  as  features  for  target  Identification.  These  patterns  for 
the  various  aircraft  can  be  seen  to  be  uniquely  different,  and  the  task 
Is  to  properly  quantify  these  differences  so  that  they  can  behave  as 
features  for  a  target  Identification  system. 

As  an  example  of  the  separability  of  the  various  aircraft  response 
curves  In  this  parameter  space.  Figure  3.9  shows  an  estimate  of  the 
smoothed  response  of  each  of  the  aircraft  plotted  on  a  single  Poincare 
graph.  Note  that  the  various  curves  are  uniquely  different.  A  small 
set  of  data  on  an  unknown  target  (from  this  group)  would  clearly  be 
Identifiable.  An  example  where  this  has  been  done  Is  shown  In  Figure 
3.10.  In  this  figure,  the  RCS  of  the  various  aircraft  was  plotted  on 
the  polar  projection  of  the  Poincare  sphere  for  the  two  orthogonal 
circular  Illumination  cases,  unly  the  starting  frequency  (S),  the 
middle  frequency  (M)  and  the  final  frequency  cases  (F)  are  Indicated. 
Note  the  wide  separation  of  the  various  parameters.  It  can  be  expected 
that  If  the  polarization  states  for  an  unknown  aircraft  were  to  be 
plotted  jn  this  curve.  It  would  be  possible  to  associate  the  unknown 
with  a  particular  aircraft  type.  (We  assume  that  the  unknown  aircraft 
Is  Included  In  the  "catalog".)  Future  research  will  focus  on  these 
concepts  In  an  attempt  to  quantify  these  results,  on  the  development  of 
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Figure  3.4.  Polar  projection  of  Poincare  sphere  representation  for 
Boeing  747  (band  limited). 
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Polar  projection  of  Poincare  sphere  representation  for 
Boeing  727  (band  limited). 
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Polar  projection  of  Poincare  sphere  representation  for 
Boeing  707  (band  limited). 
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Figure  3.8.  Polar  projection  of  Poincare  sphere  representation  for 
Concorde  (band  limited). 
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Figure  3.9.  Polar  projection  of  Poincare  sphere  representation  for  all 
'  '  five  aircraft  (band  limited). 


I 


> 


► 


» 


V 

V 

AIRCRAFT 

O 

-  CONCORDE 

V 

•  747 

o 

o 

•  727 

£ 

A 

•  707 

-0 

□ 

■  DC-IO 

NOTE:  L  s  LEFT  HAND  CIRCULAR 
R  s  RIGHT  HAND  CIRCULAR 


ILLUMINATION 
I  LLUM INATION 


Figure  3.10. 


Polar  projection  of  Poincare  sphere  representation  of  RCS 
of  five  aircraft  (selected  frequencies). 
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automatic  algorithms  to  perform  this  task  and  on  a  study  of  the 
sensitivity  of  such  techniques  to  measurement  errors  In  the  raw  data 


3.3  POLARIZATION  FEATURES  AND  SIMULATION  RESULTS 


The  Initial  Investigation  of  the  performance  of  target 
Identification  systems  using  polarization  state  data  was  based  on  the 
following  sets  of  polarization  parameters  (all  of  which  are  easily 
derived  from  the  linear-polarization  scattering  matrix  coefficients  Tyy, 

Ihh»  and  IVH^: 

Linear-polarization  coefficients 

1.  VV  Vert  leal -receive  Vertical -transmit 

2.  HH  Horl zontal -receive  Horizontal -transmit 

3.  VH  Vertical -receive  Horizontal -transmit 

Circular-polarization  coefficients 

1.  RR  Rlght-hand-clrcular-recelve  Right-hand-el rcular-transmlt 

2.  LL  Left-hand-clrcular-recelve  Left-hand-clrcular-transmlt 

3.  RL  Rlght-hand-clrcular-recelve  Left-hand-clrcular-transmlt 

Axial -Ratio  coefficients  (see  Section  3.2.1) 

1.  RAX  Rlght-hand-cl rcular-transmlt  axial-ratio 

2.  LAX  Left-hand-clrcular-transmlt  axial-ratio 

Polarization-Ratio  coefficients 

1.  VV/HH  Vertical-Vertical  relative  to  Horl zontal -Horl zontal 

2.  VH/VV  Vertical -Horl zontal  relative  to  Vertical -Vertical 

Linear  combinations  of  polarization  coefficients;  for  constants 
a,  b ,  and  c 

1.  Linear  aVV  ♦  bHH  +  cVH 

2.  Circular  aRR  ♦  bLL  +  cRL 

Concatenated  combinations  of  polarization  coefficients 


1.  Linear  VV  |  HH  j  VH 

2.  Circular  RR  |  LL  J  RL 


Using  these  polarization  parameter  sets  as  "features"  for  the 
target  Identification  algorithm,  a  series  of  simulation  experiments  were 
conducted  for  various  ranges  of  parameters.  A  description  of  these 
experiments,  along  with  a  comparison  of  the  most  Interesting  results  and 
a  discussion  of  possible  directions  for  future  research  are  given  below. 

3.3.1  Experimental  Approach 

Let  Xj^  represent  the  complex  normalized  scattering  coefficient  of 
target  k,  for  a  particular  combination  of  transmit  and  receive 
polarizations  at  transmit  frequency  fi.  (The  term  "complex  normalized 
scattering  coefficient"  refers  to  a  quantity  vrfiose  magnitude  is  the 
square-root  of  the  measured  radar  cross  section  In  square  meters  and 
whose  phase  Is  the  calibrated  phase  of  the  measured  signal.  The  radar 
cross  section  of  an  object,  and  hence  the  square  of  the  magnitude  of  the 
complex  normalized  scattering  coefficient.  Is  defined  here  In  term  of 
decibels  relative  to  one  square  meter,  l.e.,  In  DBSM.)  A  set  of  N  such 
measurements  at  discrete  frequencies  for  target  k  at  azimuth  angle  (*) 
and  elevation  angle  (v)  may  be  represented  by  the  complex  vector. 

( tp, v)  »  [XklU,v) . xkNU,v)lT  (3,6) 

A  catalog  of  test  targets  (or  prototypes)  Is  obtained  by  entering  a  set 
of  measurements  for  each  of  the  five  aircraft,  for  each  polarization 
type,  azimuth  angle,  and  elevation  angle  of  Interest  into  a  data  base 
[21.  For  this  particular  set  of  experiments,  measurements  In  the  8-58 
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MHz  scaled-frequency  band  were  available  for  each  of  the  aircraft  at  10° 
Increments  of  azimuth  angle  and  for  the  three  linear  polarization  types 
discussed  above. 

3. 3. 1.1  Noise  Model 

The  noise  model  employed  In  this  Initial  Investigation  was  the 
additive  vftlte  Gaussian  noise  model.  In  particular,  a  complex-valued 
vector  of  measurements  representing  the  measurement  of  a  radar  return, 

Y  from  an  "unknown"  target  was  simulated  by  forming  the  complex  vector 
sum 


Y  «  X  +  N  (3.7) 

of  one  of  the  set  of  cataloged  target  vectors  and  a  vector  of  complex¬ 
valued  Gaussian  deviates  N  *  [Nj,  Ng,  ....  The  "Independent" 

complex-valued  Gaussian  deviates  were  generated  by  forming  the  complex 
vector  sum,  •  A^  +  ,  vfiere  j  ■  /^T,  and  where  the  zero-mean, 

o^-varlance  Gaussian  deviates  Aj  and  B-j  were  computed  using  the  methods 
discussed  in  Abramowltz  and  Stegun  [51.  The  noise  model  is  completed  by 
the  specification  of  the  variance  (or,  average  power),  2o2  of  the 
individual  complex  Gaussian  deviates  .  For  the  present  Investigation, 
the  average  noise  power  was  specified  In  absolute  terms.  In  units  of 
DBSM.  Some  of  the  reasons  for  adopting  this  method  of  specifying  the 
average  power  of  the  additive  noise  process  are  discussed  below. 
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It  Is  common  practice  In  the  radar  literature  to  consider  the 
performance  of  a  target  detection  system  as  a  function  of  a  signal -to- 
nolse  ratio  (SNR);  thus  normalizing  the  average  power  of  the  additive 
poise  with  respect  to  the  total  power  of  the  desired  signal.  While  this 
parameter  Is  a  reasonable  "benchmark"  parameter  for  target  detection 
systems,  a  performance  evaluation  of  target  Identification  systems  based 
upon  a  specification  of  noise  power  that  Is  dependent  on  the  a-prlorl 
unknown  amplitude  of  the  measured  signal  may  be  misleading  In  some 
cases.  For  example.  It  seems  unreasonable  to  expect  an  Increased  level 
of  noise  to  be  present  during  the  Interrogation  of  a  large-cross-section 
target,  relative  to  the  level  for  a  target  with  small  radar  cross 
section. 

Other  reasons  for  this  choice  of  “normalization"  Include  the  fact 
that  signal -Independent  specification  of  noise  power  most  accurately 
represents  systems  In  cases  where  Internal  or  "front-end"  noise  Is  the 
dominant  source  of  measurement  errors.  Indeed,  the  noise  power  In  DBSM 
may  be  computed  directly  In  terms  of  the  "noise  temperature"  or  "noise 
figure"  of  the  radar  reel ever. 

In  addition,  an  absolute  reference  for  noise  power  allows  a  fair 
comparison  of  the  performance  of  systems  employing  different 
polarization  coefficients  as  features.  For  example,  the  Rlght-clrcular- 
recleve,  Right-clrcular-transmlt  (RR)  coefficient  could  be  measured 
either  directly,  using  right-circular  transmit  and  receive  antennas;  or 
Indirectly,  using  both  horizontal  and  vertical  transmit  and  receive 
antennas  and  forming  the  vector  quantity 
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Assuming  that  comparable  components  are  used  to  build  each  system,  there 
should  be  no  advantage.  In  terms  of  Immunity  to  additive  noise.  Inherent 
tn  either  design. 

The  Issue  of  fair  comparisons  Is  of  even  more  concern  when 
evaluating  the  performance  of  Identification  algorithms  based  on 
features  Involving  polarization  coefficients  that  are  nonlinear 
combinations  of  measured  quantities.  In  particular,  consider  the 
formation  of  an  axial -ratio  component  vector  from  Right -circular. 
Right-circular  (RR),  and  Right-circular  Left-circular  (RL)  measurements, 
viz. 


Clearly,  It  Is  not  appropriate  to  specify  the  average  power  of  an 
additive  noise  component  In  the  original  measurements  based  on  the 
magnitude  of  an  axial-ratio  value;  nor  would  It  be  appropriate  to 
totally  neglect  the  effects  of  “front-end"  noise,  even  If  the  measured 
signal  Is  moderately  strong. 

3, 3. 1.2  Distance  Metrics 

As  discussed  above,  one  of  the  objectives  of  this  investigation  was 
the  determination  of  the  value  of  precise  phase  Information  In  the 
Identification  process.  This  was  accomplished  by  comparing  the  results 
of  experiments  using  a  "coherent"  (complex)  distance  metric, 
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dc( Y,X)  -  IY-XI2 


(3.10) 


m 


to  those  using  a  "non-coherent"  (magnitude)  distance  metric, 

d„(Y,X)  «  1 1 Y |  -  | X 1 1 2  (3.11) 

where  |Y|  denotes  a  vector  whose  components  are  the  magnitudes  of  the 
complex-valued  components  of  the  vector  Y. 

3. 3. 1.3  Identification  Algorithm 

The  Identification  algorithm  employed  In  this  study  was  the  single 
nearest-neighbor  algorithm.  The  nearest-neighbor  algorithm  anounces 
target  Identity  1  (where  1  represents  one  of  the  targets  In  the  catalog 
set)  If 

1  =  arg  min  d(Y,Xk)  ;  k  •  1,  2,  3 . k  (3.12) 

If  the  target  identity  1  does  not  match  the  Identity  of  the  catalog 
target  used  to  simulate  the  return  of  the  unknown  target,  an  error  Is 
charged  to  the  experiment.  By  performing  many  such  experiments,  using 
independent  noise  deviates  for  each  test,  a  cumulative  sum  of  errors  may 
be  computed  for  each  test  target  (posing  as  an  "unknown"  target),  and 
for  each  level  of  average  noise  power  of  Interest. 

3.3.2  Results 

The  results  of  this  study  support  the  following  hypotheses: 

*  For  most  of  the  polarization  features  considered,  the  use  of 
precise  phase  information  produces  significant  performance 
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gains  (on  the  order  of  lOdB).  Somewhat  surprising,  however.  Is 
the  fact  that  an  Identification  error  rate  as  low  as  0.1  { 10%) 
can  be  achieved  at  a  SNR  of  about  OdB  when  phase  Information  Is 
available. 

*  Of  the  polarization  types  tested,  the  linear  co-pole  features 
(HH  and  VV)  and  features  primarily  consisting  of  linear 
combinations  of  linear  co-pole  measurements  (RCP,  LCP,  and 
RLCP)  allowed  the  most  reliable  Identification. 

*  The  tests  Indicate  that  the  overall  performance  Is  most 
error-free  at  azimuth  angles  near  90“  (broadside).  While  this 
result  is  not  surprising.  It  Is  not  totally  In  agreement  with 
the  results  of  previous  simulation  studies  that  were  based  on 
SNR  noise  power  normalizations. 

*  For  this  type  of  noise  model,  features  that  are  nonlinear 
transformations  of  (noisy)  measurements,  such  as  axial -ratios, 
are  not.  In  themselves,  reliable  features. 

It  should  be  emphasised  that  these  data  represent  only  a 
preliminary  investigation,  and  the  results  cited  herein  should  be 
considered  as  Inconclusive. 


CHAPTER  4 


COMPLEX  NATURAL  RESONANCE  EXTRACTION  FROM 
CROSS-POLARIZER  MEASURED  SCATTERING  DATA 

4.1.  INTRODUCTION 

The  complex  natural  resonances  (CNR's)  of  an  object  a^e  useful  for 
Identification  purposes.  The  success  of  this  procedure,  however, 
depends  on  careful  extraction  of  the  CNR's  as  well  as  the  Identification 
process.  In  this  chapter,  the  emphasis  Is  on  the  resonance  extraction 
process  Itself  and  we  Illustrate  the  usefulness  of  cross-polarized 
scattering  data  in  this  context.  The  availability  of  accurate, 
broadband  scattering  data  has  made  the  present  study  possible  [1],  At 
most  orientation  angles  the  cross-polarized  returns  are  at  least  ten  dB 
below  the  co-polarlzed  returns  but  the  data  are  still  useful  for  natural 
resonance  extraction.  In  fact,  a  very  dominant  specular  return  In  the 
co-polarlzed  data  usually  masks  the  resonance  effect.  This  Is  found  to 
be  the  case  for  aspects  near  broadside  for  the  aircraft.  For  both  the 
co-polarlzed  and  cross-polarized  data  the  method  of  CNR  extraction  Is 
the  rational  function  fit  method  [3],  The  target:  are  electroplated 
models  of  modern  aircraft. 


4.2.  THEORETICAL  BACKGROUND 


For  simple  objects  such  as  the  sphere  and  infinite  circular 
cylinder,  rigorous  analytical  solutions  exist  for  the  location  (In  the 
Complex  s  plane)  of  the  CNR's.  For  geometrically  complicated  targets, 
however,  the  CNR's  must  be  extracted  from  measured  scattering  data. 

Even  modest  electrical  sizes  cannot  be  treated  using  a  moment  method 
formulation  and  numerical  search  procedure. 

A  scatterer  has  an  Infinite  number  of  CNR's  and  Is  therefore  a 
system  of  Infinite  order.  The  Impulse  response  of  a  geometrically 
simple  target  or,  for  that  matter,  the  response  of  the  target  to  any 
aperiodic  Illumination,  consists  of  first  a  forced  response  as  the 
wavefront  moves  across  the  target  and  then  a  free  or  natural  response  as 
the  wavefront  moves  beyond  the  target.  As  discussed  later,  the 
separation  Is  not  distinct  for  complicated  targets.  In  the  s  (Laplace) 
domain  then  the  transfer  function  of  the  target  can  be  modelled  as  a 
residue  series  plus  an  entire  function. 

-  Cn<"> 

F(s)  ■  I  s  .  s  +  G(s,f!)  .  (4.1) 

n«l  n 

The  Cn's  are  the  residues  which  are  aspect  (n)  dependent,  the  sn's  are 
the  natural  resonances  and  G(s,n)  Is  the  entire  function.  Note  that  the 
entire  function  Is  also  dependent  on  aspect. 

Any  practical  measurement  system  Is  a  bandlimlted  system  so  only  a 
limited  number-  of  CNR's  are  significant  for  calculations  In  a  particular 
frequency  range.  The  response  to  an  Impulse  Input  Is  actually  of  the 
form 
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+  G(s,n)  ♦  N(s) 


(4.2) 
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n»l 


s  -  s. 


and 


y(t) 


n  St* 
l  Cne  n  +  g( t ,n)  +  n(t) 
n«l 


(4.3) 


-  G(s,n)  within  passband 

where  G  (s,fl)  ■  { 

0  outside  passband. 

In  a  practical  system  one  obtains  the  response  as  a  function  of  real 
frequency,  corresponding  to  Information  over  a  certain  region  of  the  ju> 
axis.  Neglecting  the  noise  for  the  moment.  If  the  unknowns  of  Equation 
(4.2),  l.e.  Cn,sn  and  G(s,n)  can  be  found  somehow,  the  equation  suggests 
that  one  has  knowledge  of  the  portion  of  the  complex  s  plane  delineated 
by  the  real  frequency  limits.  For  geometrically  complicated  targets 
there  Is  no  known  method  for  finding  the  entire  function  G(s,n). 
Therefore,  some  sort  of  approximating  function  must  be  employed  In 
searching  for  the  CNR's.  Since  the  known  data  points  are  all  on  the  ju> 
axis,  one  Is  really  performing  an  extrapolation  even  though  complete 
knowledge  of  the  entire  function  Is  not  the  goal.  For  this  reason  no 
method  of  CNR  extraction  can  claim  absolute  accuracy.  The  entire 
function  part  of  the  response,  being  unknown  and  aspect  dependent, 
causes  difficulties.  In  Prony's  method,  where  the  time-domain  response 
Is  sampled,  It  has  been  suggested  that  the  late-time  portion  of  the 
response  be  sampled  [6],  With  this  approach  the  forced  response  portion 
of  the  return  Is  avoided.  This  Is  possible  for  simple  target  shapes 
although  the  return  may  be  many  dB  down.  For  complicated  targets. 
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however,  the  response  In  the  time  domain  Is  really  a  very  complicated 
combination  of  free  response  and  forced  response.  Some  substructures  of 
the  target  are  resonating  long  before  the  Interrogating  wavefront  moves 
beyond  the  target.  The  authors  agree  with  Felsen  [71  that  subtracting 

a 

the  physical  optics  approximation  would  be  a  good  first  step  to  negate 
the  effects  of  the  *_*nt1re  function.  Unfortunately,  for  geometrically 
complicated  targets  obtaining  the  physical  optics  approximation  could  be 
a  major  effort.  We  offer  an  alternative  approach  as  outlined  below. 

The  rational  function  Is  written  as 

b0  *  DlS  ♦  ty2  ♦  ...  ♦ 

F(s)  «  Z  N 

1  +  a«.s  ♦  ♦  •••  ♦  a*. sn 

1  c  N 

rewriting  Equation  (4.5) 

N  Cn 

F(s)  -  l  - T~  +  a  +  bs 

n*l  n 

bH+l 

where  b  *  ~ —  etc, 
aN 

The  expression  a  +  bs  in  Equation  (4.6)  Is  a  simple  attempt  to  model  the 
entire  function.  It  Is  motivated  by  the  fact  that  the  Impulse  response 
of  a  number  of  simple  objects  have  Impulse  or  doublet  singularities. 
Therefore,  as  long  as  the  frequency  range  over  which  one  extracts  CNR's 
Is  not  too  largo  (the  optical  range  Is  avoided)  the  approximation 
remains  good.  The  unknowns  a  and  b  are  both  aspect  and  frequency  range 
dependent. 


»  (4.5) 


*  (4.6) 
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A  distinction  should  be  made  between  the  complex  natural  resonances 
of  the  target  and  the  zeros  of  the  denominator  of  Equation  (4.5).  The 
poles  of  Equation  (4.5)  are  not  necessarily  true  CNR‘s  as  there  are 
usually  curve  fitting  poles  which  must  be  eliminated.  For  a 
geometrically  complicated  target  the  true  CNR's  will  not  always  appear 
in  each  test  because  certain  of  the  resonances  may  be  unexcited  or  only 
weakly  excited  for  particular  orientations  of  the  target. 

Equation  (4.6)  does  not  account  for  the  different  time  delays  for 
various  components  of  the  specular  returns.  Furthermore,  these  time 
delays  are  also  aspect  dependent.  A  possible  formulation  Is  as  follows: 


F(s) 


N 

I 

n=l 


s-s. 


I 

I 

1«1 


,  -gis  -his. 

(a^e  ♦  b^se  1  ) 


(4.7) 


The  parameters  and  h^  are  real  numbers  that  are  aspect  dependent. 
Equation  (4.7)  Introduces  more  unknowns  and  Is  actually  nonlinear.  It 
may,  however,  lead  to  an  even  better  fit  of  the  measured  spectrum.  The 
Individual  time  delays  only  change  the  phase  of  the  corresponding 
specular  terms.  However,  the  time  delays  of  Equatlc  (4.7)  would  affect 
the  magnitude  of  F(s)  also.  Therefore,  the  use  of  an  algorithm  that 
requires  only  magnitude  information  will  not  offer  any  Improvement.  One 
way  to  remove  the  non-linearity  Is  to  estimate  a  priori  g^  and  hj.  This 
would  require  data  which  is  sufficiently  broadband  that  Impulse  and 
doublets  with  delay  could  be  discerned.  It  would  also  require  aspect 
and  geometrical  Information  of  the  object.  The  series  expansion  of  the 
exponential  terms  Is  not  useful  since  It  may  not  converge. 
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Details  of  the  application  of  the  rational  function  method  are 

given  In  [3].  Note  that  all  of  the  coefficients  In  Equations  (4.5)  and 

(4.6)  are  real.  Therefore  the  set  of  poles  consists  primarily  of 

complex  conjugate  pairs  with  corresponding  complex  conjugate  residues. 

% 

The  form  of  Equation  (4.6)  can  lead  to  a  different  method  of  target 
Identification  as  giver,  In  [83.  With  the  method  of  Identification  given 
In  [83,  a  and  b  play  an  Important  role  and  the  Identification  technique 
Is  aspect  dependent.  For  a  given  set  of  CNR's  and  a  given  frequency 
range,  one  first  finds  the  ranges  of  variation  of  a  and  b  for  different 
signal  to  noise  ratios.  These  parameters  as  well  as  the  CNR's  are 
stored.  Given  frequency  scattering  data  from  an  unknown  target,  the  set 
of  CNR's  are  used  to  obtain  the  a  and  b  of  Equation  (4.6)  and  are  then 
compared  to  their  corresponding  ranges  of  variation  to  make  a  decision 
on  the  Identity  of  the  unknown  target.  This  method  actually  takes 
advantage  of  the  existence  of  .he  entire  function  and  works  extremely 
well  for  high  (>15  dB)  signal  to  noise  ratios.  For  signal  to  noise 
ratios  lower  than  15  dB,  the  variation  of  a  and  b  become  too  great  for 
the  method  to  be  useful  In  Its  present  form. 

We  now  have  a  basic  understanding  of  the  theory  of  CNR  extraction 
and  are  ready  to  examine  the  results  of  CNR  extraction  using 
co-polarlzed  and  cross-polarized  scattering  data. 
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4.3.  COMPUTATION  AND  ANALYSIS 


The  specular  portion  of  the  backscatter  response  from  geometrically 
complicated  targets  Is  predominantly  In  the  co-pol arl zed  component.  For 
all  but  the  most  anomalous  target  one  can  safely  conclude  that  only  a 
weak  entire  function  component  exists  In  the  cross-polarized  return.  As 
will  be  demonstrated,  this  means  better  numerical  stability  for  the 
CNR's  extracted  from  the  cross-polarized  returns. 

The  compact  radar  cross  section  measurement  range  [11  was  used  to 
obtain  calibrated  amplitude  and  phase  scattering  data  for  the  targets. 

A  digital  filter  Is  then  used  to  numerically  time  gate  and  remove 
unwanted  clutter  from  the  time  regions  where  the  target  return  Is  known 
to  be  zero.  The  time  response  of  the  target  Is  not  altered.  Figure  4.1 
shows  a  plot  of  the  cross -polarized  frequency  response  of  an  aircraft 
target  at  broadside  Incidence  (90°).  Figure  4.2  Is  the  corresponding 
band-limited  impulse  response  of  the  aircraft.  In  Figures  4.1  and  4.2 
vertical  polarization  was  transmitted  and  horizontal  polarization  was 
received.  Figures  4.3  and  4.4  show  corresponding  frequency  and  time 
domain  plots  after  the  data  have  been  passed  through  a  digital  filter. 
The  purpose  of  the  filter  Is  to  somewhat  smooth  the  data  without 
removing  any  salient  features.  The  data  shown  In  Figure  4.3  are  then 
used  for  pole  extraction.  Equation  (4.5)  Is  used  to  fit  portions  of  the 
spectrum,  typically  applied  over  a  1.0  GHz  bandwidth.  This  Is  done 
because  of  the  extrapolation  nature  of  the  calculations.  A  wider 
bandwidth  Is  not  desirable  from  a  theoretical  viewpoint  because 
windowing  must  be  used  to  help  eliminate  pattern-fitting  poles. 
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Figure  4.4.  Cross-polarized  impulse  response  of  aircraft  at  broadside 
(90°)  after  digital  filtering. 
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Choosing  endpoints  of  the  band  at  local  minima  helps  to  Increase  the 
accuracy  of  the  solution.  Slightly  shifting  this  band,  e.g.,  using  4.0 
to  5.0  GHz  and  then  using  4.05  to  5.05  GHz  helps  In  weeding  out  curve  or 
pattern-fitting  poles. 

CNR's  have  been  extracted  for  four  Incidence  angles  0°  (nose-on), 
45°,  90°  and  180°,  using  three  polarizations  (vertical-vertical, 
horizontal -horizontal  and  vertical -horizontal ) .  The  CNR's  extracted  are 
listed  in  Table  4.1.  Only  the  oscillatory  part  of  the  CNR's  are 
listed. 

It  Is  clear  from  Table  4.1  that  no  Incidence  angle  yields  all  of 
the  extracted  poles,  but  the  cross-polarized  results  have  a  lower  rate 
of  absenteeism  than  the  co-pol arl zed  results.  Reasoning  with  regard  to 
the  specular  returns.  It  seems  clear  that  even  near  nose-on  and  tall -on 
where  the  cross-polarized  return  Is  10  dR  down  from  the  co-pol arl zed 
returns  It  Is  still  advantageous  not  to  have  to  deal  with  specular 
components.  In  the  neighborhood  near  broadside  (90°) ,  considerable 
difficulty  was  experienced  In  extracting  CNR's  from  the  co-polarlzed 
returns.  This  was  anticipated  as  It  Is  clear  that  the  entire  function 
Is  quite  large  In  this  region.  Figure  4.5,  which  shows  the  vertical - 
vertical  broadside  return  and  Figure  4.6  which  shows  the  horizontal- 
horizontal  broadside  return  are  both  relatively  flat  compared  to  the 
cross-polarized  return  In  Figure  4,3. 

Table  4.2  lists  the  average  CNR's  as  found  over  the  four  angles. 
Maximum  deviation  from  the  average  over  the  angles  are  listed  Inside 
parentheses.  Both  of  the  co-polarlzed  results  are  essentially  equally 
different  from  the  cross-polarized  returns.  Table  4.1  also  shows  that 
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TABLE  4.1 

OSCILLATORY  PARTS  OF  CNR'S  EXTRACTED  FROM  BACKSCATTER 
RESPONSES  OF  DIFFERENT  POLARIZATIONS  AND  ANGLES 


0° 

45° 

90° 

180° 

vv 

VH 

HH 

VV 

VH 

HH 

VV 

VH 

HH 

VV 

VH 

HH 

— - 

1.143 

1.017 

1 

.215 

1.134 

1.106 

1 

.259 

1.050 

1.090 

1 

.208 

1.198 

1.100 

1.408 

1.470 

1.474 

1 

.494 

1.574 

1.626 

1 

.427 

1.490 

--- 

1 

.473 

1.488 

1.439 

1.642 

1.682 

1.770 

1 

.765 

1.797 

mm  ar 

1 

.714 

1.821 

1.742 

— 

— 

— 

— - 

1.956 

1.958 

1 

.919 

— 

2.097 

— 

2.015 

2.128 

2 

.097 

2.038 

1.926 

2.330 

2.195 

2.280 

2 

.210 

2.247 

2.418 

2 

.420 

— - 

MM  M 

2 

.298 

2.341 

2.233 

2.720 

2.512 

2.592 

2 

.657 

2.656 

2.717 

— 

2.647 

2.643 

2 

.657 

— 

2.597 

MM  M 

2.905 

3.073 

2 

.965 

2.871 

— 

2 

.947 

3.065 

3 

.047 

2yQfO 

— - 

3.310 

3.287 

MM  M 

3 

.298 

3.271 

3.276 

3 

.349 

--- 

M  M  M 

— - 

3.321 

3.219 

3.586 

3.747 

3.654 

— 

3.707 

3.675 

— 

3.639 

3.781 

3 

.645 

3.572 

— 

4.026 

4.059 

4.159 

3 

.893 

4.067 

3 

.971 

4.103 

4.113 

4 

.100 

4.021 

4.145 

4.192 

4.239 

— 

— 

4.208 

--- 

— 

4 

.226 

4.241 

--- 

— 

4.718 

4.543 

4 

.627 

4.727 

4.684 

4 

.717 

4.671 

M  M  M 

4 

.596 

— 

4.510 

— 

— 

5.056 

— 

5.024 

5.136 

— 

— 

— 

— 

4.924 

4.796 

5.442 

5.505 

5. 533 

— 

— 

— 

5 

.51 

5.547 

5. 511 

— 

5.677 

5.622 

52 


TABLE  4.2 
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AVERAGE  OF  OSCILLATORY  PARTS  IN  TABLE  4.1  OVER  THE  FOUR  ANGLES 
MAXIMUM  DEVIATION  FROM  THE  AVERAGE  VALUE  IS  SHOWN  IN  PARENTHESIS 


VV 

VH 

HH 

1.227 

(.032) 

1. 

131 

(.081) 

1. 

098 

(.008) 

1.451 

(.043) 

1. 

506 

(.068) 

1. 

513 

(.113) 

1.707 

(.065) 

1. 

767 

(.085) 

1. 

’56 

(.014) 

2.008 

(.089) 

2. 

003 

(.047) 

2. 

027 

(.101) 

2.315 

(.105) 

2. 

261 

(.080) 

2. 

310 

(.108) 

2.678 

(.042) 

2. 

605 

(.093) 

2. 

637 

(.080) 

2.058 

(.061) 

2. 

928 

(.137) 

3. 

073 

3.319 

(.030) 

3. 

293 

(.028) 

3. 

248 

(.029) 

3.616 

(.030) 

3. 

666 

(.094) 

3. 

703 

(.078) 

3.998 

(.105) 

4. 

063 

(.042) 

4. 

139 

(.026) 

4.209 

(.017) 

4. 

240 

(.001) 

4. 

208 

4.647 

(.070) 

4. 

705 

(.034) 

4. 

579 

(.105) 

— - 

4, 

974 

(.050) 

4. 

996 

(.200) 

5.476 

(.034) 

5. 

576 

(.101) 

5. 

555 

(.067) 
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the  co-pol arl zed  returns  are  more  expect  dependent  In  that  more  poles  are 
missing  at  many  aspects  compared  to  the  cross -pol arl zed  returns.  As 
would  be  expected  for  a  complicated  object,  all  of  the  poles  shown  In  the 
cross-polarized  returns  are  also  present  In  one  of  the  two  co-pol arl zed 
results.  The  cross-polarized  returns  cannot  contain  any  CNR's  that  are 
not  In  any  of  the  co-polarlzed  returns. 

4.4.  USE  OF  HIGHER  FREQUENCY  DATA 

The  rational  function  method  has  been  applied  to  the  horizontal 
polarization  response  of  another  aircraft  model  at  6.0  -  12.0  GHz.  The 
CNR's  (TABLE  4.3)  appear  to  have  the  same  type  of  behavior  as  those  found 
earlier  using  1.0  -  6.0  GHz  data.  Therefore,  this  suggests  that  this 
frequency  range  Is  also  one  where  It  Is  meaningful  to  extract  CNR  and  to 
use  them  for  discrimination. 

It  Is  difficult  to  relate  these  resonances  to  structural  features  of 
the  aircraft  at  this  stage.  It  Is  possible  to  postulate  physical 
substructures  that  resonate.  However,  there  Is  no  analytical  proof. 

Also,  each  physical  substructure  would  have  more  than  one  complex  natural 
resonance  associated  with  it.  One  possible  way  to  discern  structural 
dependence  would  be  to  use  absorbing  material  on  the  various 
substructures  to  study  their  effect  on  the  CNR's, 
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4.5.  CONCLUSIONS  FROM  THE  COMPLEX  NATURAL  RESONANCE  STUDY 


The  practical  Implications  and  problems  of  CNR  extraction  have  been 
studied  using  both  co-polarlzed  and  cross-polarized  backscatterlng  data. 
The  CNR  extraction  method  and  the  required  accuracy  are  directly  related 
to  target  Identification  using  prediction-correlation.  The  CNR  locations 
must  be  quite  accurate.  The  research  reported  here  has  demonstrated  two 
main  points.  First,  the  frequency  range  which  was  used  In  previous 
studies  can  be  doubled  without  encountering  problems  with  the  CNR 
concept.  It  Is  known  that  above  some  unknown  high  frequency  limit  the 
density  of  the  CNR's  will  preclude  target  Identification  using 
prediction-correlation.  Additional  research  Is  needed  to  establish  this 
high  frequency  limit  for  practical  targets.  Second,  It  has  been 
demonstrated  that  there  are  advantages  In  terms  of  orientation 
sensitivity  and  stability  to  CNR  extraction  using  cross-polarized 
backscatter  d:ta.  The  disadvantage  Is  that  the  data  may  be  as  much  as  10 
to  15  dB  down  from  the  co-polarlzed  return.  This  point  must  be  carefully 
understood.  One  could  simply  use  cross-polarized  data  for  CNR 
extraction.  In  this  case  the  measured  scattering  data  are  controlled  and 
the  measurements  are  clutter-limited,  not  noise-limited.  The  extracted 
CNR's  could  then  be  used  In  prediction-correlation  processing  using 
co-polarlzed  and/or  cross-polarized  scattering  data.  The  Idea  of  target 
Identification  from  CNR's  using  cross-polarized  real  world  scattering 
data  needs  to  be  very  carefully  explored. 
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TABLE  4.3 

OSCILLATORY  PARTS  OF  CNR'S  EXTRACTED  FROM 
HORIZONTAL  POLARIZATION  RESPONSE  OF  ANOTHER  AIRCRAFT 
AT  6.0  -  12.0  GHZ  BETWEEN  0®  AND  90° 


15° 

30° 

45° 

60° 

75° 

90° 

6.465 

4.247 

6.245 

6.520 

6.479 

6.  bS9 

6.410 

7.406 

7.302 

7.385 

7.282 

— 

7.413 

7.342 

7.908 

8.005 

7.906 

8.000 

8.077 

— 

8.040 

8.686 

— - 

8.516 

8.556 

8.696 

— 

8.730 

8.974 

9.017 

9.141 

9.156 

9.924 

9.153 

8.918 

9.521 

9.585 

9.776 

9.550 

9.562 

9.582 

9.621 

• 

10.053 

10.117 

10.053 

•  m 

!C.Cr9 

9.951 

10.53 

10.45 

10.469 

— 

10.516 

10.475 

10.59 

10.873 

— 

11.078 

11.094 

10.95 

— 

11.072 

CHAPTER  5 


CONCLUSIONS 

A  number  of  Important  conclusions  can  be  given  at  this  time. 

First,  the  goal  of  development  of  a  target  Identification  data  base 
and  associated  data  base  management  system  was  achieved.  This  data  base 
was  specifically  configured  for  the  study  of  the  resonance-region 
broad-band  polarization  matrix  behavior  of  the  series  of  commercial 
aircraft  and  naval  ships  on  which  data  were  available  at  the 
ElectroSclence  Laboratory.  It  Is  now  possible  to  access  the  full 
polarization  scattering  matrix  (amplitude  and  phase;  co-polarlzed  and 
cross-polarized;  fully  calibrated)  of  these  ships  and  aircraft  using  a 
single  computer  file  access  system  which  permits  the  study  of  the 
behavior  of  these  various  targets  for  a  number  of  polarization  states 
and  feature  sets. 

The  polarization  matrix  data  base  was  next  used  for  the  purpose  of 
testing  a  number  of  target  Identlflcalton  concepts.  Chapter  three  has 
shown  the  details  of  the  results  for  tests  where  the  backscattered 
return  for  a  number  of  polarization  states  and  ratios  of  polarization 
states  were  used  as  features  In  the  target  Identification  algorithms. 

In  general,  It  was  shown  that  co-polarlzed  data  from  fully  coherent 
radar  systems  (1e:  systems  trfiere  the  relative  phase  of  the  returned 
signals  could  bo  measured  as  well  as  the  amplitude)  with  multi  frequency 
capability  could  perform  very  effectively.  Slgnal-to-nolse  ratios  of  0 


to  +10  OB  In  the  final  estimates  of  the  RCS  and  relative  signal  phase 
for  the  various  targets  was  shown  to  be  sufficient  for  ml sclasslfl cation 
percentages  less  than  5  percent.  (See  chapter  three.)  Of  particular 
note  was  that  the  Introduction  of  the  ability  of  the  radar  system  to 
measure  relative  received  signal  phase  would  produce  performance 
Improvements  (over  an  "amplitude-only"  system)  of  as  much  as  10  DB. 

The  use  of  complex  natural  resonances  (CNR)  In  a  radar  target 
Identification  system  was  also  studied  using  this  new  data  base  system. 
CNR  values  were  extracted  from  both  co-polarlzed  and  cross -polarized 
measured  scattering  data.  It  was  shown  that  natural  resonances  can  be 
meaningfully  utilized  using  the  6.0  to  12.0  GH2  data.  This  Indicates 
that  the  30  to  60  MHZ  band  (full  scale)  can  be  used  for  CNR  radar  target 
identification.  It  was  also  shown  that  If  the  co-polarlzed  returns  are 
dominated  by  specular  scattering  then  there  are  significant  benefits  to 
the  use  of  the  cross-pol adzed  data.  In  particular.  Increased  stability 
of  the  measured  resonances  with  respect  to  target  orientation  was 
demonstrated. 

There  are  a  number  of  Important  areas  for  future  research  which  can 
be  discussed  at  this  point.  First,  the  extension  of  these  results  Into 
higher  frequency  bands  should  be  studied.  Data  at  higher  frequencies  (a 
factor  of  two)  for  some  of  the  aircraft  exist  at  this  time  and  such  data 
should  be  studied.  In  particular,  the  application  of  these  concepts  at 
higher  frequencies  may  permit  the  use  of  such  techniques  on 
subcomponents  of  targets  rather  than  the  overall  target.  This  area  of 
research  could  become  a  major  new  area  In  Itself. 
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New  concepts  In  the  utilization  of  extracted  target  resonances 
should  also  be  explored.  In  particular,  the  application  of  alternate 
polarization  states  as  Input  to  the  resonance  extraction  techniques  may 
be  appropriate. 

New  algorithms  for  the  Identification  process  also  need  to  be 
studied.  Optimization  techniques  have  not  been  applied  as  of  yet  to  the 
radar  target  features  and  the  associated  Identification  algorithms.  The 
optimum  frequency  bands  and  polarization  states  for  Identification  of 
specific  target  types  has  not  been  derived. 
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